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Abstract 

We study the motion of a grain boundary separating two otherwise stationary 
domains of hexagonal symmetry. Starting from an order parameter equation 
appropriate for hexagonal patterns, a multiple scale analysis leads to an ana- 
lytical equation of motion for the boundary that shares many properties with 
that of a crystalline solid. We find that defect motion is generically opposed 
by a pinning force that arises from non-adiabatic corrections to the standard 
amplitude equation. The magnitude of this force depends sharply on the mis- 
orientation angle between adjacent domains so that the most easily pinned 
grain boundaries are those with an angle 4° < < 8°. Although pinning 
effects may be small, they do not vanish asymptotically near the onset of 
this subcritical bifurcation, and can be orders of magnitude larger than those 
present in smectic phases that bifurcate supercritically. 
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The microstructure of a condensed phase and the distribution of topological defects 
largely determine its mechanical and thermodynamical response, as well as the temporal 
evolution of its non-equilibrium configurations. Continuum or hydrodynamic approaches to 
phases with broken symmetry are now well understood IHQ], including the long wavelength 
description of topological defects . Nevertheless, a quantitative theory of defect motion 
(e.g., dislocation glide) remains a difficult problem because it requires short scale phenomena 
that lie beyond a hydrodynamic theory, and therefore beyond the degree of universality that 
such a description entails. 

Bifurcations to states with broken symmetry are also encountered in a variety of physical, 
chemical and biological systems driven outside of thermodynamic equilibrium. The charac- 
teristic length scales involved are much larger than atomic dimensions, therefore making the 
observation and study of defects substantially easier. Here, amplitude or phase equations 
that focus on slow modulations of a base periodic pattern play the role of the hydrodynamic 
description . It has been noted, for example, that at leading order hexagonal patterns are 
the dissipative analogues of a two dimensional, isotropic solid 0. We focus in this paper on 
grain boundary motion in hexagonal patterns, and find many qualitative and quantitative 
similarities with grain boundary motion in crystalline solids. In contrast with the latter 
case, our starting model allows a detailed analysis of the breakdown of the long wavelength 
description of defect motion through an explicit multiple scale analysis. We are then able 
to derive analytically several results that are known only qualitatively or empirically in 
crystalline solids. 

Hexagonal patterns formed in a spatially extended system (like in Langmuir monolayers 
or block copolymer melts [^) are often fairly disordered. They consist of randomly ori- 
ented grains separated by grain boundaries (arrays of dislocations in the low angle case), and 
are equivalent by symmetry to a polycrystalline solid. This nonequilibrium microstructure 
usually evolves on a very slow, fiuctuation-dependent time scale due to pinning forces to 
defect motion. Our focus here is on a coarse grained model of a hexagonal pattern, and, 
in particular, on the motion of a grain boundary separating two domains with arbitrary 
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mis-orientation. We extend to this case a recent study of defect motion and pinning near 
a supercritical bifurcation involving stripe patterns (or smectic phases) |^. We show that 
the coupling between the slow variables of an envelope description of the defect and the 
underlying (fast) periodicity of the base pattern leads to an effective periodic potential, the 
close analogue of the Peierls barrier acting on defects in a crystalline solid ||lOl- Further- 
more, we find that the magnitude of the potential barrier does not vanish as the subcritical 
bifurcation point is approached, contrary to the case of smectic phases . Hence, not only 
crystalline phases emerging from a subcritical bifurcation are harder than smectics in terms 
of the forces acting on topological defects, but also an envelope description that ignores the 
internal degrees of freedom of defects does not appear to be valid even near the bifurcation 
point. The self-generated pinning effects discussed here are expected to be a general feature 
of modulated phases, and can explain, for example, many grain boundary conformations 



that have been observed experimentally and numerically in block copolymer melts ITT . 

We study here the Swift-Hohenberg model of Rayleigh-Benard convection with an addi- 
tional quadratic term to allow the formation of hexagonal patterns 

^ = e^-l-Akl + '^'r^ + 92^'-^'- (1) 

The order parameter ilj{x,t) is related to the vertical velocity at the mid plane of the con- 
vective cell, e is the reduced Rayleigh number, and g2 can be related to deviations from 
Boussinesq behavior in the working fluid. The uniform solution ip = becomes unstable for 
e > to a periodic pattern of wavelength Aq = 2TT/kQ. Hexagonal patterns are stable for 
— I em ((72) I < e < eA/(5'2), and roll patterns for e > eA/- In this study, we chose e G [0,eAf] so 
that only hexagonal patterns are stable. 

An approximate stationary solution for a configuration containing a planar grain bound- 
ary between two uniform and symmetric domains of hexagons that have a relative mis- 
orientation angle ^ (0 < < 7r/3, see Fig. [|) can be found by assuming that ipo = 
J2n=i ^n^^^"'^ + C.C., where An{x) are slowly varying envelopes, and x is the coordinate 



normal to the boundary. As x ^ —00, Ai 2;3 Aq = {g2 + J g2 + 15e)/15 exponentially 



fast outside a boundary layer of thickness ^ [v4„(a;) = /n(a;/^)], whereas ^44^5^6 0. As 
X — 00, ^4^5 6 Aq and ^1,2,3 0. We chose \kn\ = kQ, ki 4^ = /i;o(cos(0)x=Fsin(0)y), where 
= 7r/6 — 0/2, and the other vectors kn are obtained from these by rotations of ±27r/3 as 
indicated in Figure |l|. In analogy with crystalline solids, a small angle grain boundary is well 
described as an array of penta-hepta defects (the dislocations cores of a hexagonal pattern), 
separated on average by a distance of the order Xo/6- However, since the projections of ^1,2,3 
on the X axis are usually not commensurate, the patterns observed along the boundary are 
not periodic. 

We next focus on grain boundary motion and extend our earlier results for stripe patterns 
0]. The amplitude equation formalism has already been used to study defect dynamics in 
hexagonal patterns |T2|. Amplitude equations are obtained from a standard multiscale 
analysis of Eq. ([1|). For example, the equation for Ai follows from the solvability condition 

^ hm ^[L(^o) + 92^Pl - ^!\e-'''-'' = 0, (2) 
where L is a linear operator that follows from Eq. (|l]) |jl3| and Ix a length of 0(Ao) to be 



specified later. When both e and (72 0, the length scale of variation of the A^s {i.e. the 
grain boundary thickness ^) is much larger than Aq, and only non oscillatory terms contribute 
to the integral (0). The solvability condition leads to the usual amplitude equations that 
show that defects are either immobile or move with constant velocity 0. 

If, on the other hand, e and (72 are small but finite, then the amplitudes are not strictly 
constant within a period l^. Any oscillatory term in of wavevector K parallel to the 
normal a;-axis does not integrate to zero, and the equation for the slowly varying amplitudes 
cannot be decoupled from the phase of the defect (yielding "non-adiabatic corrections" 
14-|T6t3). Although these contributions are small (non-analytic in both e and (72), they may 



have dramatic effects. In the case of Eq. (|T]), non- vanishing terms can arise from the cubic 
nonlinearity, and are proportional to AfA4^, AiA^, AiA^A^, AiA-^A^, A1A2AQ and A1A2AQ. 
An oscillatory term proportional to cos{Kx') in Eq. (|^) contributes to order exp{—\K^\) 
to the law of grain boundary of motion (see Eqs. (^ and (0) below). This contribution 



is simplest when both e ^ 1 and (72 ^ 1 so that ^ Aq and a single mode with lowest 
K dominates. In this limit, the slowest non-adiabatic corrections are those proportional to 
A1A2AQ and A1A2AQ, which have K = \k2 + k^l = 2A;osin(^/2) (see Fig. p. The solvability 
condition (H) now reads 

dAi OF 



dt dAi 

Jxo Ix 

where = Ao/[2 sin(0/2)]. The first two terms in Eq. represent the standard amplitude 
equation with the Lyapunov functional 



df 
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- 2^2(^1^2^3 + A^A^Aq + A1A2A3 + A^A^Aq) 

+ n E \An\ + 3 E l^nl \A-m\ 
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with Dn = d/dxn the derivative along Equations similar to (j^) can be derived for the 
remaining amplitudes An (not shown). The last term in the r.h.s. of Eq. (^) is new and 
represents the dominant non-adiabatic correction in the limits e -C 1 and (72 ^ 1- 

In order to derive an equation of motion for the grain boundary from the equations for the 
Ani we first denote hy an{x) {1 < n < 6) the leading order amplitudes of the stationary grain 
boundary, solutions of the one-dimensional coupled Ginzburg-Landau equations dF/ddn = 
0. We then look for solutions of the form An{x,t) = a„(x — Xgb(t)), where Xgb{t) is the now 
time-dependent position of the grain boundary 0. We find, 

F) Xgb = -Phex sin[2/coa;gbSin(6'/2)], (5) 

where D = (ix X]n=i(^a;'^n)^ is a friction coefficient, and 

c/x COS [2A;o sin(26')x + /3] I 3 a2d:i:{alaQ) 
+ a69^(aga2)] + 12 [ai9^ (010206) + 03^2; (030205) 
+ 04(9^(040205) + 05^2; (05O2O6)]} (6) 



is the (dimensionless) amphtude of a pinning force of wavelength Ao/[2 sin(^/2)]. This pe- 
riodic force explicitly arises from the coupling between the slowly varying wave amplitudes 
and the periodicity of the base state in the integral term of Eq. (0). Equations (^) and 
(^) show that a planar grain boundary initially located at an arbitrary position relaxes 
toward the nearest minimum of the periodic pinning potential. The stationary and sta- 
ble positions of the boundary are thus discrete, separated from each other by a distance 
Aa^gb = Ao/[2 sin(^/2)] > Aq- The two wavevectors with the smallest projection on the 
grain boundary normal {k2 and k^, see Figure ^ set the wavevector {k2 + k^) of the peri- 
odic pinning potential. The usual amplitude equation formahsm would instead predict that 
Phex = 0, and Xgb is arbitrary and decoupled from the phase of the pattern. Figure |I] shows 
three successive stable locations of the the grain boundary obtained by numerically solving 
Eq. (|I|). The values of Axgb determined numerically agree very well with the analytic result. 

We expect equations of the form of Eq. (|^) to be generic at a bifurcation, and not 
limited to hexagon-hexagon grain boundaries. In particular, we anticipate that non-adiabatic 
effects are important in block copolymer melts, and can explain the conformations of planar 
interfaces observed in these systems . 



Remarkably, the present result for a pattern of hexagonal symmetry is analogous to the 



Peierls force acting on dislocations in crystalline solids . Peierls calculated the energy of a 
single dislocation by summing over the interactions between atoms within two neighboring 
layers, their displacements given by continuous elasticity as a first approximation. The 
energy of the dislocation oscillates as a function of its position so that it can glide only if a 
force of finite amplitude acts on it (the Peierls' force). Here, we find that a similar force acts 
on assemblies of dislocations organized in arrays. Like in the simpler geometry studied by 
Peierls, the amplitude of the pinning force decays exponentially with the spatial thickness 
of the defect ^. From Eq. (H), we find 

with c* and a* dimensionless constants of order unity that can a priori depend on the 



mis-orientation 6. At this point, we note an important distinction between defect motion 
in hexagonal patterns that emerge at a subcritical bifurcation, and in stripe patterns that 
bifurcate supercritically. The supercritical case was discussed in ref. for a 90° boundary 
separating two domains of stripes given by Eq. ([l|) with g2 = 0. The corresponding pinning 
force acting on the boundary, Pstripe, satisfies a relation similar to Eq. (0), with ^ ~ 
Hence Pstripe ~ exp(— l/y^) 0, as the control parameter e ^ 0. In the hexagonal phase, 
however, ^(e = 0) ~ 15Ao/(8-\/67r5f2) is finite. Therefore Phex, although small, does not 
vanishes when e — (nor when e — — |em(5'2)|)- 

Figure ^ shows typical variations of the pinning force as a function of e for different 
values of = ^r/Q). The curves are given by Eq. (^, where the amplitudes a„ have 

been obtained by numerically integrating the system of equations dF/ddn = 0. For a value 
of g2 as small as 0.3, phex is many orders of magnitude larger than Pstripe- Therefore, non- 
adiabatic effects in hexagonal, "crystalline-like" , patterns are difficult to avoid. Defects need 
to overcome much higher activation barriers, and will be much less easily un-pinned either 
upon the application of external stresses or random noise (that would be represented by 
additional terms in the r.h.s. of Eq. (^). 

Finally we study the dependence of the pinning force with grain boundary mis-orientation 
9. Figure ^ displays Phox(6') for different values of the parameters g2 and e. We find that 
the least mobile grain boundaries, i.e. those for which the amplitude of the pinning force is 
maximal, have a low angle 6m, typically such that 4° < 6m < 8°. Both 6m and the overall 
shape of Phex{6) depend weakly on g2 and e. Figure |^ also shows the pinning force vs. 6 
obtained from a direct numerical solution of Eq. (0) with g2 = 0.3 and e = 0.05. Here phcx 
is estimated by fitting grain boundary relaxation trajectories to Eq. (|^), with the friction 
coefficient D assumed to be given by the analytic result. The numerical results compare 
reasonably well with the theory given the numerical difficulties in tracking the relaxation 
of the grain boundary in a finite sized system. Although the maximum value is lower and 
the curve flatter in the numerical case, the right order of magnitude is obtained, as well 
as a similar range for 6m- The discrepancy can also be attributed in part to non-adiabatic 
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corrections of higher order to Eq. (|^), that may not be neghgible at g2 ~ 0.3. 

Our present analysis can quahtatively explain the properties of polycrystalline, partially 
ordered configurations that are typically observed at long times in many pattern forming 
systems with this symmetry Although the evolution defined by Eq. (|l]) is driven by 

the minimization of a Lyapunov functional, we expect that grain boundaries (and other 
topological defects) will become pinned at long times as the driving force for microstructure 
coarsening decreases. Therefore, the system will eventually reach metastable, glassy-like 
configurations that can only order by slow activated processes (presumably logarithmic in 
time, as already observed in ^ with a random forcing term added to Eq. (|l|)). We addi- 
tionally note that defects observed in cold metals are essentially low angle grain boundaries 
in the range 5° and 10° [|T^, a value that compares very well with the values obtained by 
the present theory. 

In summary, we have analyzed the motion of grain boundaries in hexagonal patterns from 
an order parameter equation and extended the standard Ginzburg-Landau equation for the 
slowly varying amplitude to incorporate non-analytic corrections. Like in crystalline phases, 
defect motion is opposed by short range forces with periodicity and amplitude that strongly 
depend on the mis-orientation angle between domains. Although small, these pinning forces 
can not be neglected asymptotically at long times in a coarsening system, even near onset, 
and are orders of magnitude higher than those produced in patterns of smectic symmetry. 

We are indebted to Frangois Drolet for useful discussions. This research has been sup- 
ported by the U.S. Department of Energy, contract No. DE-FG05-95ER14566. 
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FIGURES 




FIG. 1. Order parameter ip in gray scale. Only a portion of a square grid of 512^ points is 
shown with spacing Ax = \/3/2 and Aq = 8Ax. Three locations are shown in which a planar 
grain boundary is stationary. The two hexagonal domains are defined by the sets {/^i, A;2, ^3} and 
{/c4, /js, fee} respectively, with a mis-orientation angle 9 = 7r/9. 
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FIG. 2. Values of the pinning force as given by Eq. (^) as a function of e for various values of 
92 (6* = 7r/9 in each case). The result corresponding to a stripe pattern was given in ref. |^], and 
is shown as a comparison. 
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FIG. 3. Angular dependence of the pinning force for various values of §2 and e obtained ana- 
lytically with Eq. (^), and numerically for one of the parameter sets. 
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